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Summary 

This  report  presents  the  development  of  a  likelihood  ratio  test  for  synchronous  and 
balanced  CDMA  signals.  The  hypothesis  associated  with  a  binary  test  is  the  code  length 
and  number  of  user.  The  approach  incorporates  a  statistical  dependency  of  the 
spreading  matrix  elements  when  computing  the  expectation  over  the  spreading  matrix 
elements.  The  dependency  is  a  discrete  Gaussian  probability  function  that  depends  on 
a  modified  Total  Square  Correlation.  The  resulting  average  likelihood  is  a  product  of 
terms  that  depends  on  the  matrix  that  achieves  the  lowest  Total  Square  Correlation.  If 
the  set  of  signals  is  constrained  to  Hadamard  matrices  using  Sylvester's  construction, 
the  likelihood  function  can  be  simplified  further,  allowing  the  classification  over  code 
lengths  in  the  range  of22  to  213  and  possibly  higher. 

Keywords:  DS/CDMA  signals,  classification,  balanced  CDMA  load,  synchronous  CDMA, 
decision  theory,  average  likelihood  ratio  test,  Total  Square  Correlation 
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1.  Introduction 


Modulation  classification  is  an  important  step  in  the  characterization  of  signals  and 
the  recovery  of  the  information  content.  In  a  non-cooperative  scenario,  an 
eavesdropper  attempts  to  gain  information  on  how  the  signal  has  been 
synchronized  and  modulated  before  trying  to  recover  the  information  content.  A 
second  scenario  could  be  the  case  of  a  smart  receiver  with  software  defined  radio 
that  reconfigures  itself  for  various  modulation  types,  managing  the  allocation  of 
users  and  adjusting  the  code  length  for  efficient  utilization  of  the  bandwidth. 

The  problem  of  signal  classification  has  been  widely  studied  for  MPSK  and  QAM 
signals  using  a  decision  theoretic  approach.  Surprisingly  there  is  a  lack  of  research 
on  the  area  of  classifying  CDMA  signals  using  average  likelihood  methods.  [1] 

Traditional  approaches  for  signal  classification  include  heuristics,  feature-based 
classification  and  decision  theory.  The  development  of  decision  theoretic 
approaches  must  result,  by  virtue  of  the  theory,  in  an  optimal  classifier  that 
guarantees  the  lowest  probability  of  error  in  classification.  If  there  were  such 
classifier  that  could  perform  better,  then  we  would  either  omit  an  important  piece  of 
information  in  our  assumptions  or  have  discovered  a  novel  approach  that  performs 
better  than  the  current  state  of  the  art  in  the  classification  theory.  The  idea  of 
developing  average  likelihood  methods  in  CDMA  sounds  attractive  because  it  could 
also  result  in  better  estimation  methods.  In  a  future  research,  we  would  consider 
the  ambitious  goal  of  estimating  the  CDMA  spreading  matrix  or  even  estimating  the 
information  contents  (multi-user  detection)  without  knowledge  of  the  spreading 
matrix. 

The  objective  of  this  research  is  to  develop  a  decision  theoretic  approach  for 
modulation  classification  of  CDMA  signals.  This  approach  will  be  based  on  the 
classical  average  likelihood  ratio  test  (ALRT).  We  begin  this  classification  problem 
by  finding  the  size  of  the  spreading  matrix  that  generated  the  DS-CDMA  signal.  As 
the  number  of  unknown  variables  grows,  the  averaging  process  becomes  an 
extremely  complex  task.  In  the  multiuser  detection,  a  closely  related  problem, 
authors  use  Singular  Value  Decomposition  (SVD)  to  estimate  the  unknown  variables 
and  do  estimation  based  on  a  generalized  likelihood  approach.  [2]  These 
approaches  are  suboptimal  by  nature. 

The  benefits  of  the  theory  sound  promising;  however,  there  is  a  high  risk  due  to  the 
difficulty  of  averaging  over  a  large  number  of  parameters  and  obtaining  an  exact 
analytical  solution.  During  the  first  year  of  the  research,  the  risk  was  mitigated  by 
concentrating  in  the  development  an  empirical  likelihood  function  for  CDMA.  The 
understanding  was  possible  by  using  symbol  algebra  software  and  attempting  to 
solve  simple  cases  of  CDMA.  This  knowledge  was  collected  in  a  previous  Interim 
Technical  Report  AFRL-RI-RS-TR-2013-058.  Most  of  the  previous  work  was 
superseded  by  the  research  done  during  the  second  year. 
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In  the  second  part  of  the  research,  an  effort  was  made  to  improve  and  provide  a 
mathematical  foundation  for  deriving  an  empirical  likelihood  function  found  during 
the  first  year.  This  includes  extending  the  likelihood  function  to  higher  dimensions 
of  the  spreading  matrix.  A  cluster  of  48  computers  was  used  calculate  the  likelihood 
to  code  length  dimensions  5  and  6.  The  work  required  translating  an  empirically 
derived  Matlab  code  (see  Appendix  G)  to  a  Multiple  Programming  Interface  in  C++. 
Unfortunately  the  resulting  approach  was  computationally  overwhelming  for  code 
lengths  above  6.  However,  the  code  length  cases  of  5  and  6  provided  more  evidence 
of  the  dependency  of  the  likelihood  on  matrices  that  achieve  the  lowest  Total  Square 
Correlation.  At  this  point,  it  was  clear  that  the  average  likelihood  function  for  CDMA 
is  a  sum  of  products  of  two  main  terms,  that  is,  some  coefficients  that  depended  on 
the  Total  Square  Correlation  and  a  product  of  hyperbolic  cosine  functions.  This 
product  came  from  an  averaging  of  exponentials  as  shown  in  Appendix  E.  But  using 
this  product  would  require  partitioning  the  set  of  all  possible  combination  of  the 
spreading  matrix  in  carefully  selected  subsets  as  shown  in  Appendix  C.  Another 
interesting  observation  was  made.  Most  of  the  coefficients  of  the  likelihood  function 
are  negligible  for  computing  the  likelihood  as  shown  in  Figures  3  and  4.  This 
property  reduced  greatly  the  complexity  of  the  average  likelihood.  The  final  task 
was  to  provide  experimental  evidence  that  validates  our  development.  A 
classification  using  Receiver  Operating  Characteristic  (ROC]  curves  revealed  that 
the  approach  taken  in  this  research  was  able  to  produce  desirable  results. 

This  report  presents  a  mathematical  development  of  the  CDMA  average  likelihood 
function.  It  is  divided  in  the  following  parts:  a  review  of  existing  techniques,  a 
description  of  the  model,  the  problem  and  assumptions,  the  formulation  of  the 
likelihood  function  and  its  simplification.  The  metric  of  performance  used  in  the 
simulations  is  the  ROC  curves  which  provide  a  visual  characterization  of  the  new 
classification  rules.  The  results  appear  to  support  the  validity  of  the  newly 
generated  classification  rules. 

2.  Theoretical  Background 

The  classification  of  DS/CDMA  signals  should  not  be  confused  with  the  problem  of 
multiuser  detection.  The  multiuser  detection  deals  with  the  estimation  of 
information  symbols  in  a  CDMA  system.  This  can  be  done  in  the  absence  of  the  code 
that  generated  the  signal.  The  different  approaches  try  to  gain  information  of  the 
spreading  matrix  by  using  Singular  Value  Decomposition  (SVD)  methods.  A 
common  metric  is  the  bit  error  rate  (BER).  On  the  other  hand,  the  goal  of 
modulation  classification  is  the  selection  of  a  modulation  scheme  over  the  others. 
The  classifier  can  be  characterized  using  Receiver  Operating  Characteristic  (ROC) 
curves  or  by  generating  confusion  matrices. 

As  an  estimation  problem,  the  number  of  users  or  information  symbols  is  obtained 
by  maximizing  a  conditional  likelihood  function.  Common  assumptions  in  this 
problem  are  the  estimation  under  Additive  White  Gaussian  Noise  (AGWN)  and  the 
full  synchronization  of  the  CDMA  signal.  This  includes  frame,  chip  and  phase 
synchronized  system  [3].  Maximizing  the  likelihood  function  is  prohibitively 
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complex  due  to  the  large  amount  of  unknown  parameters  that  are  required  to  be 
estimated. 

There  is  a  lack  of  likelihood  methods  for  modulation  classification.  A  survey  on 
classification  methods  discusses  only  traditional  modulation  schemes  such  as  MPSK 
and  QPSK  methods  [4].  A  more  recent  paper  highlights  the  need  for  modulation 
classification  [5],  but  none  of  the  methods  directly  addresses  the  classification  of 
CDMA  signals.  A  feature  based  classification  technique  is  shown  presented  in  [6]. 
This  approach  classifies  the  signal  by  its  code  length.  It  uses  cyclostationary  features 
and  a  neural  network  for  classification.  The  method  requires  training  and  there  is 
no  guarantee  of  achieving  an  optimal  performance. 

3.  Synchronous,  Balanced  CDMA  Average  Likelihood 
3.1  Model 

The  goal  is  to  classify  CDMA  signals  by  their  code  length  L  and  number  of  users  U  in 
the  presence  of  Additive  White  Gaussian  Noise  (AWGN).  A  CDMA  signals  of  size 
Lx  1  (see  Eq.  1)  is  generated  by  a  spreading  matrix  C  of  size  L  x  U,  and  the 
information  symbol  vector  b  of  size  U  x  1.  For  simplicity,  we  assume  that  the  signal 
the  energy  per  symbol  E  of  each  user  is  constant,  i.e.,  a  balanced  load.  Developing 
the  case  of  unbalanced  load  is  possible,  but  not  considered  in  this  report  due  to  time 
constraints. 


x  =  JIT/L  ■  Cb 

(1) 

The  spreading  matrix  and  the  information  vector  are  binary  antipodal.  The 
spreading  matrix  is  a  block  matrix  composed  of  column  vectors  c*  j  with  low  cross¬ 
correlation.  A  metric  that  characterizes  the  spreading  matrix  is  known  as  the  Total 
Square  Correlation  (TSC).  A  modified  TSC  is  defined  in  this  paper  by  [7] 

u- 1 u- 1 

t (C)  A  ^  XlMl- 

i= 0  7=0 

j*i 

(2) 

An  alternative  definition  of  [Eq.  2]  for  an  antipodal  matrix  can  be  expressed  in  terms 
of  the  Frobenius  Norm  as: 


t(C)  A  \\CT C\\2f  -  L  ■  U . 

(3] 

The  noise  vector  n  has  zero  mean,  and  covariance  matrix  £  =  N0/2  ■  /,  where  /  is 
an  L  x  L  identity  matrix.  The  receiver  has  no  knowledge  of  the  spreading  code,  the 
code  length  or  the  number  of  users.  The  channel  is  modeled  by: 
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y  =  x  +  n. 

(4) 

This  model  makes  the  following  assumptions:  chip  synchronization  and  knowledge 
of  the  beginning  of  the  sequence.  For  simplicity,  our  approach  will  use  similar 
assumptions  to  those  used  in  multiuser  detection,  i.e.,  chip  synchronization  and 
frame  synchronization.  The  assumptions  represent  a  scenario  where  the 
eavesdropper  has  information  of  the  preamble  and  knows  the  beginning  of  the 
transmission. 

The  conditional  likelihood  function  is  given  by: 

A(y|?f  =  {L,U}.C.b)  =  I 

(5) 

The  likelihood  is  conditioned  on  the  hypothesis  Of  =  {L,  U),  the  spreading  matrix  C, 
and  the  information  vector  b  according  to  (5).  The  likelihood  has  made  the 
assumption  of  the  code  length  L  and  frame  period.  The  assumption  must  be  tested 
against  other  possible  assumptions  in  order  to  make  a  determination  of  the  code 
length. 

The  assumption  of  frame  synchronization  is  not  a  critical  restriction  when 
constructing  an  average  likelihood  function.  Additional  averaging  will  be  sufficient 
to  remove  the  assumption.  Once  it  is  removed,  it  will  represent  a  more  realistic 
scenario  where  the  eavesdropper  does  not  know  the  beginning  of  the  sequence. 

Under  the  assumption  of  no  frame  synchronization,  our  conditional  likelihood  is 
given  by: 


X{y\0f  =  {L,  {/},  C.b.e) 


1  (ye-x)HX  1{ye-x) 

ONo)t/2  6 


(6) 


where  ye  represents  the  received  vector  that  is  delayed  e  chips. 


3.2  Average  Likelihood  Function 

The  construction  of  our  likelihood  function  will  be  expressed  as  a  function  of  the 
energy  per  chip  EC  =  E/L  and  noise  power  N0.  Equation  (7)  introduces  new 
variables  for  convenience.  These  are  the  energy-to-noise  ratio  per  chip  y,  the 
unnormalized  signal  s,  and  the  correlator  output  r  respectively: 


y 


E 

N0L' 


s  =  Cb,  r  = 


y 

V^T 2' 


(7) 


The  conditional  likelihood  function  of  (5)  becomes: 
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(8) 


KW.C.i)  = 


(See  Appendix  A)  Unlike  the  MPSK  likelihood  [1],  all  terms  in  (8)  are  kept  because 
of  their  dependency  on  the  hypothesis. 

Up  to  this  point,  the  discussion  has  not  added  any  new  development.  The  new 
contribution  is  contained  in  the  following  premise.  The  average  of  the  conditional 
likelihood  is  an  interesting  case  of  dependency  between  the  elements  of  the 
spreading  matrix.  It  is  unlikely  that  a  random  selection  of  matrix  elements  would 
result  in  a  low  correlation  matrix.  Therefore,  we  must  selectively  choose  the 
spreading  matrices  based  on  some  code-design  metric.  This  approach  constructs  a 
Gaussian-like  discrete  probability  distribution  based  on  our  definition  of  the  TSC. 
The  probability  of  choosing  a  suitable  spreading  matrix  with  low-correlation  is 
given  by: 


P(C)  = 


e-/?T(C) 


(9) 


For  a  perfectly  uncorrelated  code  such  as  Walsh  codes,  the  TSC  metric  is  zero  and 
the  probability  achieves  a  maximum  value.  This  probability  has  multiple  maxima. 
Any  sign  inversion  or  permutation  of  rows  and  column  of  matrix  C  must  have  the 
same  TSC  value  and  thus  multiple  maxima  are  expected.  (See  Appendix  C) 

The  parameter  /?  is  called  precision  and  it  is  half  the  inverse  of  the  variance  of  this 
particular  distribution.  [8]  Consider  p(0-C)  being  the  probability  of  a  hypothesis  J-C. 
The  detection  of  K  would  depend  on  p(fK)  integrated  from  a  threshold  p  to  infinity. 
The  resulting  probability  of  detection  as  a  function  of  p  is  a  sigmoid  function  called 
error  function.  If  /?  is  the  precision  of  p(J-C )  and  this  variance  is  allowed  going  to 
infinity,  then  the  slope  of  the  sigmoid  becomes  steeper  to  a  point  when  the 
probability  of  detection  is  either  zero  or  one.  This  adjustment  in  precision  results  in 
a  hard  decision  scheme.  Later  on,  we  will  investigate  the  case  when  the  precision  of 
this  probability  approaches  to  infinity.  Averaging  over  the  code  coefficients 
eliminates  the  dependency  of  the  likelihood  on  the  spreading  matrix  as  follows: 

A(r\H,b)  =  (g£)  c„  P(C) 

(10) 

Finally,  the  dependency  on  b  is  removed  by  averaging  over  uniformly  distributed 
information  symbols  as  given  by: 


Approved  for  Public  Release;  Distribution  Unlimited 
7 


A(r\K)  = 


/  e~^  \1/2 

IwJ 


Z  ^  Z P(C)' 


?A/2y(r,s)-y(s,s) 


bj 


All  Ci j 


(u) 


3.3  Simplification 

Averaging  process  for  hypothesis  K  =  {2,2}  has  8  unknown  variables  and  produces 
2 8  =  64  exponential  terms.  In  general,  the  number  of  possible  combinations 
involved  in  the  averaging  process  increases  exponentially: 


2  L-U+U 


(12) 


For  a  4x4  spreading  code,  the  number  of  averaging  computations  is:  242+4  =  220. 

The  likelihood  expression  in  {11}  is  difficult  to  solve  analytically.  One  can  represent 
the  different  combinations  of  information  symbols  bj  and  elements  q ;  using  a 
decimal  number  representation  that  are  converted  into  binary  digits  by  (13}.  (See 
Appendix  B) 


2lu-  1 

A(r\K)=  ^  P(Z(m))-A(f|Jf,Z(m),l) 

771  =  0 


eA/2r2f=olRe{ri>2y=o1Zi,j(Tn)-yEiio12y=o1<z*,i,z»,j>-/?-rsC(Z(m)) 

Z(m)  =  {. ztJ}LxU 

zi,i  =  Cijbj 


ZU0 n)  =  (_i)14™’<bS42) 

(13} 

However,  additional  simplification  comes  from  exploiting  the  invariance  properties 
of  the  TSC  and  energy  of  the  CDMA  signal  to  rows,  columns  and  sign  inversion 
transformations  of  the  spreading  matrix.  Using  invariance  transformations  can  be 
useful  in  simplifying  generalized  likelihood  ratio.  [8]  The  space  of  possibilities  of 
the  Z  matrix  can  be  partitioned  in  two  subsets  and  <SgT  defined  below. 

Definition  1:  The  subset  is  the  subset  of  matrices  Y  =  {yq  }  and  Z  =  {zjj} 
that  meet  the  following  condition: 
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U- 1 

U-l 

X - 1 

i? 

II 

y  yu 

= 

X  Zij 

VI 

o  1 

7=0 

Eq.  14 

Lx  f/  LX  U 

Definition  2:  S^  T  is  the  subset  of  matrices  Y  =  { y^- }  and  Z  =  { zt  j )  such  that: 


t/— 1 

x — 1 

U-l 

x — 1 

Clj  — 

y 

= 

y 

VI 

o  1 

i-i 

TSC(Y )  =  TSC(Z )  =  r 


(15) 


The  set<Sa;T  is  a  subset  of<Sa  by  definition.  The  TST  and  energy  are  invariant  under 
row/column  permutations  and  row  sign  inversions. 


The  subsets  S$  are  disjoint  sets  of  all  possible  combinations  of  the  matrix  C.  We  can 
partition  the  each  set  S g  in  disjoint  subsets  S^T  as  illustrated  in  Figures  1  and  2. 


c 

^ a 

Co 

L. _ 

Figure  1.  Partition  Using  Sets  Sa 


Figure  2.  Partition  Using  Subset  S g 


Proposition  4:  Let  Z  be  a  matrix  subset  !PC  a  matrix  permutation  with  binary 
elements  {+1  or  0},  and  Ts  a  diagonal  matrix  that  consists  of  ±1  entries.  Then  the 
transformation: 


Z'  =PS-Z-TC 


(16) 


is  also  contained  in  the  subset  S$  T.  (See  proof  in  Appendix  C.) 

We  reformulate  (11)  in  terms  of  summation  terms  using  Definitions  1  and  2  in  (14) 
and  (15)  as  follows: 
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A(r|Jf)  =  ^  ^  ^  ^  P(Z)  ■  A(r|Jf,Z,  1). 


a  r  ZE  ,%T  Tin  J>s 


(17) 


It  can  be  shown  (see  Appendices  D  and  E)  that  the  average  likelihood  becomes: 

L— 1 

A(r|Jf)  =  ^  aL:UiS(p)  J~[  e~Yl>J=^ cosh(Jly  ■  Refo]  ■  at) 

All  a  i= 0 

2l 

2  Count(ai=0)^All  T^tZE  S-£T  & 

<*L,u,a(P)  = - zl - : - — ■ 

^iAll  a~fcount{af=^)^AllT^Z€  S-^T^  ^ 

(18) 

The  likelihood  function  can  be  expressed  as  a  product  of  cosh  functions.  The 
argument  of  each  cosh  function  depends  on  the  sum  of  the  rows  of  matrix  Z.  The 
determination  of  all  the  subsets  S^r  introduces  is  extremely  difficult  because  it 
requires  fully  characterizing  2L  U  spreading  matrices  in  terms  of  <SgT.  There  is  no 
analytical  formula  available  for  such  characterization;  however,  studying  the 
behavior  of  aLU^ (/?)  can  lead  us  to  find  an  alternative  simplification. 

3.4  Exact  Solution  of  a  2x2  Spreading  Matrix 

The  exact  likelihood  expression  can  be  computed  using  Mathematica.  A  desktop 
computer  can  calculate  the  likelihood  for  code  lengths  2  <  L  <  4  as  shown  in 
Appendix  F.  Alternatively,  we  can  compute  ccL  U  ^(p)  using  the  Matlab  code  provided 
in  Appendix  G. 

The  likelihood  for  a  2x2  spreading  matrix  hypothesis  is  given  by: 

/  e-{f,f)  \1/2 

M?\H)  =  L  )2  j  («2,2,[o,o]t  cosh(0  ■  g0)  cosh(0  ■  g±)  + 

a2,2,[2,o]Te~4y  cosh(2  ■  g q)  cosh(0  ■  gx)  + 
a2,2,[o,2]Te~4y  cosh(0  ■  g o)  cosh(2  ■  gx)  + 
a2,2,[2,2]Te~8y  cosh(2  ■  go)  cosh(2  ■  g x)) 
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where 


a2,2,[0,0]T  ~  a2,2,[2,2]T 
a2,2,[2,Qi]T  =  a2,2,[0,2]r 


1  Z 

21  +  z 
1  1 
21  +  z 


z=e 


z=e 


9i  =  J2 y  ■  Rein}. 


[19] 


The  behavior  ofai(/^(/?)  will  be  characterized  in  terms  of  infinite  precision,  i.e., 
lim^oo  aL  U  z  (/?)  .  The  convergence  to  zero  provides  an  important  simplification  of 
(20).  Several  ccL  U  ^(p)  parameters  have  been  calculated  for  code  lengths  up  to  4  as 
shown  in  Appendix  H.  The  parameters  that  converge  to  a  non-zero  value  have  been 
identified  with  an  asterisk  in  the  Appendix. 

Definition  3.  The  vector  a*  =  [a0,  is  a  feature  vector  if  C*  =  {ci,j}LxU  is  a 

spreading  matrix  that  achieves  the  minimum  TSC  and: 

\u-i  | 

f=o  I 


Property  5.  The  limit  of  ctLU^(P)  as  /?  approaches  infinity  depends  on  the  feature 
vectors  as  follows: 


lim  aLyUiS  (/?) 

/?->oo 


ro 

2L  I  c.  I 

<  2  C°unt(ai=0)  |  ^  CL,Tmin\ 

y  —  i  c . 

cl  '2^Cou7it{cip=-0')  |  d’^TTiin 


for  a  ^  a* 


for  a  =  a* 


(20) 


The  proof  is  shown  in  Appendix  J.  For  ft  oo,  the  aL  U$  parameter  is  just  the  ratio  of 
the  cardinality  of  our  previously  defined  sets  <SgT  and  <Sg  multiplied  by  a  constant. 
Increasing  the  precision  is  analogous  to  deriving  a  K-Means  classifier  from  the 
Gaussian  Mixture  Models.  [8]  Figures  3  and  4  show  plots  of  cc3  2ia(P)  in  logarithmic 
scales.  It  can  be  appreciated  that  a3:2,a(P)  converges  very  rapidly  for  /?  >  1.  Similar 
observations  were  made  for  code  length  cases  2  <  L  <  6.  (See  Appendix  I)  For  code 
length  of  higher  length,  we  would  have  to  make  the  assumption  that  these 
observations  are  also  valid;  otherwise,  the  ocL  Uiz(P)  coefficients  would  need  to  be 
calculated.  This  is  already  a  computationally  intractable  problem. 
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For  antipodal  matrices,  we  provide  the  following  empirical  rule  for  computing  the 
parameters  aL  U  First,  we  find  any  matrix  C*  with  minimum  TSC.  Then  we 

multiply  the  matrix  by  all  possible  2U  choices  of  the  vector  b  and  take  the  absolute 
value  of  each  element  of  the  resulting  vector.  The  result  is  a  set  of  feature 
vectors  a*.  The  limiting  value  of  aL  U  ^  is  the  proportion  of  a  given  vector  a  in  the  set 
of2y  possibilities.  (See  Appendix  J) 


lim  aLUa*  (fi)  = 

/?->oo 


count(abs(C*b )  =  a*) 
2" 


for  all  possible  b. 


(21] 


-0.9 

-L.0 


3  i'  2  +  6  Z  +  Z 1 


4  (6  4-  9  z  +  zJ 


*3 


Lo£l^_y]  I.i3|| 
-14  C~ 


3  <2  +  z) 

4  (6  +  9  z  +  z1) 


IP\ 


Figure  3.  Non  Vanishing  Coefficients  for  3x3  Spreading  Matrix 


*13  lip] 

Lo£tdrj!  [£11  3  3  1011 


4  (6  +  9  z  +  z1) 


Figure  4.  Vanishing  Coefficients  for  3x3  Spreading  Matrix 

In  order  to  prove  this  assertion  (21),  we  need  to  demonstrate  that  two  antipodal 
matrices  C1  and  C2  with  same  minimum  TST  and  different  sum  of  rows  are  always 
related  by  rows  and  column  permutations  C1  =  P1C2T2.  This  property  is  referred 
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as  a  Hadamard  equivalence,  but  unfortunately  it  does  not  always  hold.  [10] 
Therefore,  it  would  be  necessary  to  find  all  possible  C*  that  are  not  equivalent  or 
constrain  our  set  to  a  specific  construction.  In  our  case,  we  will  constrain  our 
problem  to  the  Sylvester’s  Construction  of  the  Hadamard  matrix.  The  coding  theory 
does  not  provide  yet  an  answer  for  finding  an  expression  of  the  lowest  TSC  matrices 
of  an  arbitrary  size. 


3.5  Average  Likelihood  Function  Example 

In  this  section  we  are  going  to  summarize  the  steps  needed  for  constructing  the 
likelihood  function  given  the  hypothesis  K  =  {4,4}. 


Step  1:  Find  a  matrix  with  minimum  TSC 


This  step  is  easy  for  Hadamard  matrices  with  L  =  U  =  2k. 


C*  = 


1 

1 

1 

-1 


1 

-1 

1 

-1 


Step  2:  Generate  the  feature  vectors  using  process  described  in  [Eq.  21]: 
dl  =  [4, 0,0, Of,  aLtUiSl  =  2/24  =  1/8 

=  [2.2,2,2]T,  y  g*  =  8/24  =1/2 


Step  3:  Form  the  product  of  hyperbolic  cosine  functions  for  the  different 
permutations  of  the  feature  vectors. 

Permutations  of  a\\ 

[4, 0,0, Of,  [0,4,0, Of,  [0,0, 4, Of,  [0,0,0,4f 
For  each  one  aL  U  =  2/24  =  1/8 
Permutations  of  a^: 

[2,2,2,2f 

with  ocL  Ufi*  =  8/24  =  1/2 


Product  of  cosh  functions: 


M? =  {4,4})  =  y^nNoyj  e-42ycosh(4  ■  1  Rein})  + 

1  3  \ 

-  e_22ycosh(2  ■  J2 y  ■  f?e{rj)  J 


Approved  for  Public  Release;  Distribution  Unlimited 

13 


As  the  code  length  and  number  of  user  grow,  the  average  likelihood  becomes 
complicated,  but  the  complexity  is  much  less  than  the  original  formulation  found  in 
(11).  The  220  unknown  variables  have  been  totally  removed  from  the  original 
expression.  In  the  next  section,  the  expression  will  be  simplified  by  adding  a 
constraint  to  the  problem. 

3.6  ROC  Using  Simplified  Likelihood 

Figures  5  through  14  show  ROC  curves  for  several  simple  cases  of  spreading  codes. 
The  length  of  the  signal  under  test  is  expressed  as  the  total  number  of  chip.  This 
value  is  based  on  selecting  5!  ■  2  =  240  symbols  of  the  hypothesis  with  the  highest 
code  length.  If  the  code  length  is  2  (see  Figure  5),  length  of  the  signal  is  480  chips. 
The  ROC  curves  are  plotted  for  different  SNR  values.  All  of  the  simulations 
correspond  to  /?  -»  oo.  Cases  3~C  =  { L ,  1}  are  trivial  ones.  Their  likelihood  ratio 
reduces  to  a  constant  when  compared  to  a  BPSK  signal  and  it  is  impossible  to 
classify  them  correctly  as  shown  in  Figure  6. 


Figure  5.  ROC  curves  for  ALRT  Test 

Jf1  =  {2,2}vs.3f0  =  {l,l} 


Figure  6.  ROC  curves  for  ALRT  Test 
3f1  =  {2,l}vs.Jf0  =  {l,l} 


Figure  7.  ROC  curves  for  ALRT  Test 

U  i  =  (3, 3}  vs.  JC  o  =  (1, 1} 


Figure  8.  ROC  curves  for  ALRT  Test 

=  (3, 3}  vs.  K o  =  (2, 2} 
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Figure  9.  ROC  curves  for  ALRT  Test 
Jt,  =  {3, 3}  vs.  Jt0  =  {3, 2} 


HOC C uv* U sm g  'jin pH ic -Aw 
HD.  H^3i;-3S-3-3  W.  h  I .  BHSK-SiHM 


Figure  11.  ROC  curves  for  ALRT  Test 

=  {4, 4}  vs.  K0  =  {3,3} 


Figure  13.  ROC  curves  for  ALRT  Test 
X,  =  {4,4}vs.2f0  =  {2,2} 


Figure  10.  ROC  curves  for  ALRT  Test 

=  {3, 3}  vs.  It  o  =  {2, 2} 


Figure  12.  ROC  curves  for  ALRT  Test 
It,  =  {4, 4}  vs.  Jt0  =  {3,2} 


Figure  14.  ROC  curves  for  ALRT  Test 
It,  —  {4,4}vs.  Jf0  =  {1, 1} 


The  orange  curves  are  extreme  cases  where  the  chip  SNR  is  so  low  that  the  classifier 
starts  breaking.  As  more  noise  is  added  to  the  signals,  the  ROC  curves  invert  their 
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performance.  A  noisy  CDMA  signal  under  hypothesis  null  hypothesis  Jf0  will  look 
like  K±.  This  might  be  a  result  of  numerical  precision  errors. 

The  developed  classifier  is  a  hard  decision  scheme.  We  have  defined  the  joint 
probability  of  the  spreading  code  in  terms  of  P(C),  but  we  can  also  define  the 
probability  in  terms  of  the  TSC  as  P(r|C;  /?).  As  it  is  expected  from  hard  decision 
classifiers,  the  probabilities  becomes  discrete  instead  of  continuous:  P(r  = 
T min\C)P )  ~  1  and  P(t  gt  Tmin\C]p)  «  0. 

We  might  expect  that  codes  with  a  slightly  deviation  from  the  minimum  TSC 
spreading  matrices  would  be  detected  and  classified  correctly.  The  main 
parameters  affecting  the  performance  of  the  classifier  are  the  parameter  /?  and  the 
feature  vector  a*.  If  we  refer  to  Figure  4,  we  can  notice  that  for  a  small  precision 
such  as  /?  >  1,  the  coefficients  «L,t/,a(/?)  are  already  small  enough  to  consider  valid 
the  approximation  of  Equation  (20).  So,  the  parameter  that  mostly  affects  the 
performance  of  the  classifier  must  be  the  parameter  vector.  If  two  codes  has  similar 
parameter  vector,  then  the  classifier  rule  should  also  work.  This  is  only  an 
assessment.  Future  work  must  consider  the  scenario  of  small  deviations  from  the 
minimum  TSC  codes. 

3.7  Classification  of  CDMA  Signals  based  on  Sylvester's  Construction 

We  consider  the  simplified  hypothesis  K  =  { L ,  L }  using  the  Sylvester’s  Construction 
of  the  Hadamard  matrix. 

Definition  4:  A  Hadamard  matrix  H  is  a  L  x  L  matrix  of  elements  ±1  that  satisfies 
the  following  condition  [11]: 


HtH  =  L  ■  II 


(22) 


/L  is  an  L  x  L  identity  matrix. 


Definition  5:  The  Sylvester’s  Construction  of  the  Hadamard  matrix  is  given  by  the 
following  recursive  formula  using  the  Kronecker  product  [10]: 


h2 


+  1 
-1 


] 


H2k  —  H2®H2k-i 


(23) 

Hadamard  matrices  have  been  a  subject  of  study  of  mathematicians  for  years.  There 
exists  several  constructions  known  at  the  present  time;  however,  there  is  no 
theorem  that  proves  the  existence  for  any  arbitrary  code  length  [8]. 
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It  is  possible  to  generate  distinct  Hadamard  matrices  using  the  Kronecker  product  of 
two  Hadamard  matrices. 

Theorem  2:  The  Total  Squared  Correlation  of  a  matrix  generated  from  a  Kronecker 
product  is  given  by:  [12] 

TSC(Ha®Hb )  =  TSC(Ha)TSC(Hb )  +  LaUaTSC(Hb )  +  LbUbTSC(Ha ) 

[24] 


Corollary  2-1:  The  Kronecker  product  of  two  Hadamard  matrices  is  another 
Hadamard  matrix. 

For  a  Hadamard  matrix  the  TSC  is  zero  based  on  Equation  [3].  The  expression  in 
[24]  is  zero  under  the  Kronecker  product  of  two  Hadamard  matrices. 

Corollary  2-2:  The  Total  Square  Correlation  of  a  Hadamard  matrix  is  zero. 

Unfortunately  the  Kronecker  product  of  two  minimum  TSC  matrices  does  not 
necessarily  produce  a  matrix  with  minimum  TSC  unless  both  matrices  are 
Hadamard. 

The  feature  vector  corresponding  to  a  Sylvester’s  Construction  is  given  by: 


=  [fo  =L.Si  =  0 . h-i  =  Of 


[25] 


as  well  as  the  L  possible  row  permutations  of  the  feature  vector  f*.  This  is  a  fact 
known  from  the  Sylvester’s  Construction  Hadamard  Matrix. 

By  restricting  our  set  0  to  S |*  0  we  eliminate  other  constructions  of  Hadamard 
matrices  by  having  setting  our  ratio  to  be  one: 


2l 

sr,o 

2Count(^i=0) 

2l 

^?*,0 

2Count^i=0) 

=  1 


[26] 


The  average  likelihood  function  of  such  scheme  reduces  to  the  following  expression: 
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— (r,r)  \  2 

T^Ny)  Z  J- 1  e~r^il2  cosh(A/f2y  ■  ti  ■  /?e{rj}). 

K  '  Ml  i=o 

Permutations 

of  f  *  elements 

[2  7) 


For  the  case  L  =  1,  the  expression  reduces  to  the  average  likelihood  function  of  a 
BPSK  signal  found  in  [1]. 

The  construction  of  the  likelihood  ratio  between  hypotheses  K2L  =  H2l  and 
J-CL  =  Hl  provides  additional  simplification.  (See  Appendix  K)  For  every 
vector  r2Lxl  in  K2L  there  are  two  vectors  rLxl  infCL.  The  two  vectors  are 
statistically  independent,  so  the  likelihood  in  the  denominator  is  multiplied.  The 
terms  e~('r’r'>  containing  the  energy  of  the  vectors  r2Lxl  and  rLxl  are  cancelled.  The 
final  form  is  given  by: 

_ _ Xfoip  cosh (^2/  •  2 L  •  fie{rfo}) _ 

ACrIJfi)  cosh  (fly  ■  L  ■  Re{rk })  ■  Z^l+i  coshCT^y  ■  L  ■  i?e{rfe}) 


(28) 


Using  the  double  angle  formulas  for  cosh  functions  helps  reducing  the  amount  of 
calculations  by  reusing  the  terms  in  the  denominator. 

A(r|Jf2[)  ZgL0(2 dg  -  1) 

2llrJ  A(?|Jft)  a,o  d*  •  Sl.1+1 4 

=  cosh(1/2y  ■  L  ■  i?e(rfc}) 


(29) 


The  product  of  cosh  functions  can  cause  numerical  overflow  in  the  computations. 
The  likelihood  ratio  shown  in  the  previous  equation  offers  better  numerical 
stability.  Nevertheless,  this  simplification  implies  that  full  Hadamard  CDMA  signals 
have  easily  implementable  rules  and  can  be  easily  detected. 

4.  Results 

The  ratio  test  in  (28)  allows  computing  the  likelihood  ratio  for  matrices  of  higher 
code  length  with  low  complexity.  The  Figures  15  through  26  show  the  likelihood 
ratio  calculated  for  two  sets.  The  length  of  the  signal  under  test  is  expressed  as  the 
total  number  of  chip.  This  value  is  based  on  selecting  5!  ■  2  =  240  symbols  of  the 
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hypothesis  with  the  highest  code  length.  If  the  code  length  is  2  (see  Figure  15), 
length  of  the  signal  is  480  chips. 


Figure  15.  ROC  Curves  Figure  16.  ROC  Curves 

Hadamard  4  versus  2  Hadamard  8  versus  4 


Figure  17.  ROC  Curves  Figure  18.  ROC  Curves 

Hadamard  16  versus  8  Hadamard  32  versus  16 


Figure  19.  ROC  Curves 
Hadamard  64  versus  32 


Figure  20.  ROC  Curves 
Hadamard  128  versus  64 
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Figure  21.  ROC  Curves 
Hadamard  256  versus  128 


ROC  Curve  Using  Simplification 
HO: Had  512  vs.H1:Had  1024 


llsinn  7457fi0rhin<i;5pmnlp  1000  4am nips 


Figure  23.  ROC  Curves 
Hadamard  1024  versus  512 


ftOCQra  Uirig  'iinpUicalioo 
HO  Hud3«aw.HI.H5d40S6 


Figure  25.  ROC  Curves 
Hadamard  4096  versus  2048 


Figure  22.  ROC  Curves 
Hadamard  512  versus  256 


ROC  Curve  Using  Simplification 
H0:Had  1024  vs.H1:Had  2048 


Figure  24.  ROC  Curves 
Hadamard  2048  versus  1024 


Add  Cl m  Uwg 

HO.  Had  ^0%  w. HI  Hofl  61 12 


Figure  26.  ROC  Curves 
Hadamard  8192  versus  4096 


The  chip-level  SNR  appears  to  be  extremely  low.  First,  we  note  that  the  chip-level 
SNR  is  a  quantity  defined  for  convenience  and  can  be  expressed  in  terms  of  a  more 
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meaningful  quantity:  the  symbol-level  SNR.  The  symbol-level  SNR  is  the  energy  to 
noise  power  density  ratio  of  one  symbol: 

E 


dB(ys )  =  dB(yc )  +  3.01  ■  k. 


m 

The  symbol  SNR  equals  the  chip  SNR  plus  3  dB  times  the  code  length  L  of  the 
Hadamard  matrix.  The  symbol  SNR  should  not  be  confused  with  the  total  SNR  of  the 
CDMA  signal. 

The  low  SNR  can  be  explained  by  understanding  the  process  of  generating  a  CDMA 
signal.  When  a  full  load  CDMA  signal  is  constructed,  there  is  a  probability  psp;fce  for 
finding  one  single  spike  of  magnitude  L  and  energy  L2  within  a  frame  of  L  —  1  zeros. 
The  signal  classifier  can  be  interpreted  as  a  detector  of  spikes  with  probability  pspjfce 
equals  to  Equation  (21). 


5.  Conclusions 

1.  This  research  demonstrated  that  is  possible  to  construct  a  classifier  for  the 
detection  of  full-loaded  CDMA  signals. 

2.  The  simplest  CDMA  classifier  is  a  hard  decision  classifier  that  only  considers 
the  matrix  that  achieves  the  lowest  Total  Squared  Correlation  for  a  given 
code  length. 

3.  Further  simplification  of  the  likelihood  equations  can  be  made  if  a  particular 
construction  of  the  spreading  code  is  assumed.  Such  is  the  case  of  Sylvester’s 
construction  of  Hadamard  matrices. 


This  approach  defines  the  probability  of  a  low  correlation  matrix  and  uses  this 
definition  for  averaging  over  all  possible  spreading  matrices.  The  resulting 
likelihood  is  a  sum  of  the  product  of  two  terms  as  shown  in  (18).  The  first  term 
depends  on  the  Total  Square  Correlation  and  a  parameter/?  and  the  second  term 
depends  on  the  feature  vectors  a*.  The  first  term  acts  like  a  filter  that  dampers  the 
effect  of  many  of  the  products  of  cosh  functions.  After  dropping  the  insignificant 
terms,  the  expression  terms  depend  exclusively  on  the  spreading  matrix  that 
achieves  the  minimum  Total  Square  Correlation. 
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The  detection  of  Hadamard  matrices  generated  using  the  Sylvester’s  Construction 
adds  more  simplification  to  the  likelihood  function.  The  sum  of  the  product  of  cosh 
functions  reduces  to  only  one  summation  as  shown  in  (29).  Although  the  expression 
is  in  a  simple  form,  numerical  overflow  is  a  concern  for  any  implementation  of  the 
likelihood  function  using  large  code  lengths. 


6.  Future  Work 

The  case  of  balanced  CDMA  signal,  i.e.,  equal  energy  signals  provides  a  limited 
application  of  the  theory.  The  next  step  in  this  research  will  be  the  assumption  of 
different  energy  levels  for  each  signature  vector.  Also,  the  case  of  QPSK  CDMA  will 
be  investigated.  Preliminary  research  has  shown  that  it  is  possible  to  develop 
similar  equations  in  terms  of  the  product  of  cosh  functions  and  feature  vectors. 


Approved  for  Public  Release;  Distribution  Unlimited 
22 


7.  Bibliography 


[1]  C.  Huang  and  A.  Polydoros,  "Likelihod  Methods  for  MPSK  Modulation 
Classification,"  IEEE  Transactions  on  Communications,  1995. 

[2]  Y.  Yao  and  V.  Poor,  "Eavesdropping  in  the  Synchronous  CDMA  Channel:  An  EM- 
Based  Approach,"  in  IEEE  Transactions  on  Signal  Processing,  2001. 

[3]  Y.  Yao  and  V.  Poor,  "EM-Based  Blind  Demodulation  of  Synchronous  CDMA," 
IEEE  6th  Sym.  On  Spread  Spectrum  Tech,  and  Applications,  2000. 

[4]  0.  Dobre,  A.  Abdi  and  W.  Su,  "Survey  of  Automatic  Modulation  Classification 
Techniques:  Classical  Approaches  and  New  Trends,"  IET  Communications,  2007. 

[5]  J.  Xu,  W.  Su  and  M.  Zhou,  "Likelihood  Ratio  Approaches  to  Automatic 
Modulation  Classification,"  2011. 

[6]  S.  Khamy,  A.  Helw  and  A.  Mahdy,  "Classification  of  Multirate  CDMA  Signals 
using  Compressed  Cyclostationary  Features,"  Progress  in  Electromagnetics 
Research  Proceedings,  2012. 

[7]  P.  Cotae,  "Distributive  Iterative  Method  for  Minimizing  Generalized  Total 
Square  Correlation  of  CDMA  Systems,"  IEEE  ITW,  2003. 

[8]  C.  Bishop,  Patter  Recognition  and  Machine  Learning,  Springer,  2006. 

[9]  P.  Schreier,  L.  Scharf  and  A.  Hanssen,  "A  Generalized  Likelihood  Ratio  Test  for 
Impropriety  of  Complex  Signals,"  in  IEEE  Signal  Processing  Letters,  2006. 

[10]  "Encyclopeadia  of  Design  Theory,"  School  of  Mathematics  Sciences,  Queen 
Mary,  University  of  London,  [Online].  Available: 
designtheory.org/library/encyc.  [Accessed  16  September  2013]. 

[11]  E.  Assmus  and  J.  Key,  "Hadamrad  Matrices  and  their  Designs,"  Transactions  of 
the  American  Mathematical  Society,  vol.  330,  no.  1, 1992. 

[12]  P.  Lancaster  and  M.  Tesmenetsky,  Theory  of  Matrices  with  Applications, 
Academic  Press,  1985. 

[13]  J.  Hamkins,  "Low  Complexity  Approximation  to  Maximum  Likelihood  MPSK 
Modulation  Classification,"  in  Jet  Propulsion  Laboratory,  National  Aeronautics 
and  Space  Administration,  2004. 


Approved  for  Public  Release;  Distribution  Unlimited 

23 


[14]  F.  Hameed,  0.  Dobre  and  D.  Popescu,  "On  Likelihood-Based  Approach  to 
Modulation  Classification,"  IEEE  Transactions  on  Wireless  Communications,  vol. 
8,  no.  12,  2009. 

[15]  A.  Vega  Irizarry,  "Optimal  Methods  for  Classification  of  Digitally  Modulated 
Signals,"  AFRL/DTIC,  2013. 


Approved  for  Public  Release;  Distribution  Unlimited 
24 


APPENDIX  A.  Likelihood  Form 


We  define  the  following  variables: 


Signal  to  noise  ratio  per 
chip 

E 

Yc  ~  n0l’ 

Unnormalized  signal 

s  =  C  -b 

CDMA  transmit  signal 

II 

Output  of  the  correlator 

.  y 

Then,  we  proceed  to  substitute  them  in  the  likelihood  function: 


1  {y-x)HI.-1(S-x) 

A(y \K  =  { L,U},C,b )  =  ,  ..  ,,<0e 


A(y\K  =  {L,UlC,b)  = 


(7TAi0)L/2 


OW/2 


2N0/2 


1  -^rHr+2^Re{yfN^/2  rHCb]-j-(CYb)HCb 


No 


A(y|^  =  {L,l/},C,h)  = 


(jzNq)1^2 


The  final  result  is: 


e  2 


__  1/2 

A(y|Jf  =  {1, !/},  C,  b )  =  eVWM?")™ 


-+y2yfte{fH}C£>-y(C£>)  Cf> 


-y(Cb)H  Cb 
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APPENDIX  B.  Average  Likelihood  Function 


Consider  the  base  2  and  decimal  representation  of  b  and  C  combinations  as  follows: 

nio  =  ■■■  >B1,B0)2 

mio  =  (Cl-u-ij  ■■■ ,  Ci,  C0)2. 

The  multiple  summations  over  binary  symbols  bj  and  q j  can  be  implemented  with  a 
summation  of  the  decimal  indexes  n  and  m  respectively: 

bj(n)  = 

c,j  (m)  = 

2LU  —  1  2U  —  1 

A(r\M)  =  ^  ^  P(C(m))  ■  A(f|^f,C(m),b(n)) 

771  —  0  71  —  0 

A  new  variable  zi j  equivalent  to  the  product  bj  and  q;-  eliminates  the  need  of  one 
summation. 


Zij(.P)  =  ctJ  ■  bj 


(_1)1+7nod(l^vJ’2) 


The  likelihood  reduces  to: 


where: 


2LU —1 

A(r \K)  =  ^  P(Z(m))-A(r\H,Z(m),l) 

M  =  0 


A(f\0-C,Z(m),l) 


ey[Zyl%=o 


Z(jYl)  iz*,0>z*,l>  ->z*,U-l\ 


and 


1  =  [1 . If- 


Note  that  the  metric  is  invariant  to  this  transformation  because  |b[|2  =  1. 
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U- 1 U- 1 

r(C)  ^  ^  Xl<^^)|2 

i=0  7=0 
j±i 


t(Z,  b ) 


2 


f/—l  t/-l 

r(Z)  =  ^  ^J<z*)i,z*j}|2  =  t(C) 

i=0  7=0 
i*i 
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APPENDIX  C. 


Invariance  to  Permutations 


LxU  LXU 

Definition  1:  The  subset  is  the  subset  of  matrices  Y  =  { ytj }  and  Z  =  {zjj} 

that  meet  the  following  condition: 


U- 1 

x — 1 

U- 1 

x — 1 

di  — 

y. 

= 

V.  z« 

VI 

o  1 

l-l 

Eq.  31 

LxU  LxU 

Definition  2:  S$  T  is  the  subset  of  matrices  Y  =  {yij}  and  Z  =  { zitj }  such  that: 


U- 1 

x — 1 

U- 1 

x — 1 

II 

y  y>j 

= 

X  z« 

1 

II 

o  1 

1  O 

II 

TSC(Y )  =  TSC(Z )  =  r 


Eq.  32 


Proposition  1:  The  TSC  is  invariant  to  permutations  and  sign-inversions  of  rows 
and  columns. 

Proof:  A  permutation  matrix  T  has  the  following  property: 

T-Tt  =  1 

Using  our  definition  of  a  TSC  for  a  binary  antipodal  matrix  [Eq.  4], 

TSC(?1C?2)  =  WTl^Tf^C^W}  -L-U 
=  W^TlT^TfWl-L-U 
=  \\CtC\\2f  -L-U 
=  TSC{C ) 

Proposition  2:  The  TSC  is  invariant  to  transposition. 

Proof:  The  Frobenius  norm  is  invariant  to  transposition.  The  term  L  ■  U  does  not 
change  after  transposing  the  matrix  C. 
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As  an  important  note,  it  is  interesting  to  mention  that  Hadamard  matrices  are 
known  to  be  invariant  to  permutations  of  rows  and  columns,  sign  inversions  of  rows 
and  columns  and  transposition.  [10] 

Proposition  3:  The  energy  of  a  signal  is  invariant  under  the  following 
transformation  using  permutation  matrices: 


b'  =  Txb 
C'  =  T2CTf 

Proof:  This  is  an  easy  proof  using  permutation  matrices.  The  energy  of  a  CDMA  per 
symbol  signal  is  given  by: 


xTx  =  ybTCTCb 


x'Tx'  =  yb'T  C'T  C'b' 

x'Tx'  =  yfa  bf(3>2  C?l)T(?2  C?l )  fa  b) 
x'Tx'  =  yb'T?l?1CTPf?2C?lT1b 
x!Tx!  =  ybTCTCb 
x'Tx'  =  x Tx 

Proposition  4:  Let  Z  be  a  matrix  subset  <%T,  Tc  a  matrix  permutation  with  binary 
elements  {+1  or  0),  and  Ts  a  diagonal  matrix  that  consists  of  ±1  entries.  Then  the 
transformation: 

Z'  =TS-Z-TC 

Eq.  33 


is  also  contained  in  the  subset  S^  r. 

Proof:  It  is  assumption  thatZ  £  £gT.with: 


CLi  = 


,  and 


TSC(Y )  =  TSC(Z )  =  r. 
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Also,  it  is  claimed  that  Z'  e  S$T.  under  Z'  =TS-  Z  •  Tc.  Therefore, 


CLi 

U- 1 

=  \ps-z-tc-i\  = 

U- 1 U- 1 

^  ^  "  zi,j  (Pc)  j,k  '  1 

j=0 

k= 0  7=0 

U-l 

f/-l 

CLi  = 

= 

Z Zij 

7=0 

7=0 

The  expression  (Tc)j,k  does  not  change  the  sign  of  the  column  z*  ;-.  It  only  sorts  the 
columns.  Also, 


TSC(Z)  =  TSC(Z')  =  t 


by  Proposition  1.  Therefore,  under  the  transformation  Z'  =  Ts  ■  Z  ■  Pc,  the 
transformation  is  also  part  of  the  subset:  Z'  G  S$  T. 
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APPENDIX  D.  Simplification  of  the 
Likelihood  Function 


The  next  step  in  the  simplification  of  the  likelihood  requires  partitioning  the  matrix 
space  Z  in  subsets  of  equal  TSC  and  energy.  Averaging  over  the  likelihood  will  be 
replaced  with  an  average  over  matrices  belonging  to  all  subsets  and  all  the 
possible  sign  inversions. 

=  Z  Z  Z 

Alla  z'eSnAll?s 

An  element  of  a  that  equals  zero  introduces  double  counting  of  spreading  matrices 
in  the  subset  <Sg.  To  correct  this,  we  must  divide  by  2  raised  to  the  sum  of  all  zero 
elements  in  a. 


The  entire  summation  can  be  calculated  if  the  whole  set  is  known  for  every 
possible  vector  a.  The  summation  over  the  sign  inversion  matrices  Ts  can  be 
simplified  as  follows: 


^  P{Z)-X{r\M,Z,T) 

All  ?s 


z 

All  pj  j=+ 1 


2i= o  ^eiriiPi,i^=  o  zi,j 


Using  Theorem  1  of  Appendix  E,  results  in 


z 

All  Piit=±  1 


o  Re{ri)Pi,i  HjL o  zi,j  — 


z  ■■■  z 

P  o,o=±1  Pl-i,l-i=±1 


eZi=o  Pi,i(V2rReird2y=  o1  Hj) 


L— 1  U- 1 

=  2L  Y [  cosh (y/2r  ■  n  ^  Zij) 
i= 0  7=0 


L—l 


=  2L  Y [  cosh (VZy  ■  r£  ^ 

i=0  7=0 


1/-1 


) 


L— 1 

=  2L  J~[  cosh(72y  ■  rt  ■  a*) 

i= 0 


So, 

L-l 

^  P(Z)  ■  A(f|Jf,  Z,l)  =  ]~|  coshC^/Zy  ■  Re{rJ  ■  a*) 

v4ZZ:Ps  i=0 
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L—l 

=  2Le~P'T(z)  J^e_yla‘l2cosh(A/27  ■  Refo}  ■  a*) 

i= 0 

The  sum  over  can  be  expressed  in  the  sum  of  all  possible  TSC  values  r  and  all  the 
subsets  Sa.T- 

A(r\K)  =  I  2Sum(ai=0)  III  P(Z)  ■  A(r\H,Z,T) 

Alla  All  r  Z  G  S-£T  All  Ts 


£  P{Z)-A(r\H,Z,T) 


=  ~I 


All  Tc 


V 


L—l 


U- 1 


2Sum(ai=0)  /  , 

Alla  AllrZeS 


2le  J- 1  e  Y^i=o  IzmI  cosh(i/2y  ■  Refo}  ^  z t  j) 

eSa,r  i=0  7=0 


=iy _ i _ y  y 

W  Z_i  2  Sum(ai=0) 


/lllfl 


All  T  Z  G  <£7; 


L—l 

2Le~P'T  J~ |  e-ya^ cosh(A/2y  ■  Refo}  ■  a*) 

i= 0 


-z(n 


?  7^'=o  ai  cosh(A/2y  ■  Re{rJ  ■ 


e 


-/?T(Z) 


w  = 


y  2L  y  y  e-w^ 

/  .  2 5um(aj=°)  /_! 

a  AUtZES-£t 


The  likelihood  of  a  CDMA  signal  can  be  expressed  as  a  sum  of  a  product.  The  first 
term  depends  only  on  a  and  can  be  expressed  as  a  product  of  cosh  functions.  The 
second  term  depends  on  a  and  r.  The  variable  matrix  Z  can  be  replaced  with  the 
spreading  matrix  C.  This  equation  is  consistent  with  the  exact  likelihood  calculated 
for  simple  CDMA  cases  as  shown  in  Appendix  F. 
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APPENDIX  E.  Discrete  Exponential  Averages 


Theorem  1:  Let  x  and  u  be  vectors  in  RLxl.  The  average  of  the  exponent  of  the  dot 
product  between  x  and  u  over  the  vector  elements  x*  =  ±at  is  given  by: 


L—l 

\  ^  ^  =  ]~[cosh(Uj  -at). 

x0=±a0  x,.-i=±a,.- 1  i= 0 


XL-i=±aL- 


Proof:  This  property  can  be  proven  by  induction.  Consider  vectors  x!  and  u'  be 
vectors  in  RL+lxl.  The  dot  product  of  vectors  x'  and  u'  can  be  constructed  by 
multiplying  by  the  exponential  term  eUL'XL  on  both  sides  of  the  equation: 


*o=±ao 


1 

'2 


I 

XL-l=±dL-l 


L—l 

eUTX+UL-XL  _  euL-xL  n  cosh(iq  ■  at ) 

i= 0 


Then,  averaging  over  the  variable  xL  results  in: 


L—l 

I  -  I  eu  X+ULXL  I  eUL-XL  n  cosh (Ui  ■  at ) 

x0=±a0  xj-±aT,  xT-±aT.  i= 0 


xl=±cll 


xl=±cll 


L—l 

=  cosh(eUL'XL)  J” |  cosh(i^  ■  at ) 

i= 0 
L 

=  ]~|  cosh(iq  ■  af) 

i= 0 

Therefore,  the  formula  applies  to  vectors  in  RL+1. 


Corollary  1-1: 


iz 

x0=±a0 


\  I 

xl-i=±o-l-i 


L—l 

cosh(urx)  =  PJ  cosh(tq  ■  a*) 

i= 0 


Proof:  The  proof  is  derived  from  Theorem  1.  It  requires  decomposing  the  cosh 
function  in  a  sum  of  exponentials. 


iz 

x0=±a0 


I  Z 

^L-l=±«L-l 


cosh(urx) 


1  l  euT*  +  e  uTx 

2  Zj  "'2  Zj  2 

x0=±a0  *l-i=±ol-i 


5Co  =  ±aO  xL-l  =  ±aL-l  X0=±a0  xL-l  =  ±aL-l 
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Corollary  1-2: 


1  L_1  1  L_1 
=  —  J- 1  COSh(iq  ■  CLj)  +  —  J- 1  COSh(tq  ■  CLi ) 

i=0  i= 0 


L— 1 


=  J~[  cosh(Uj  ■  at) 

i= 0 


\  2  sinh(^) = 


x0=±a0  xL-1=±aL- 


Proof:  The  proof  is  derived  from  Theorem  1.  It  requires  decomposing  the  cosh 
function  in  a  sum  of  exponentials. 


—  ^  ^  eUTx=~2  '"2  cosh(urx)  +  sinh(urx) 

x0=±a0  xL-1=±aL-1  x0=±a0  xL-1=±aL-1 


L—l 

=  ]~ [  cosh(Uf  ■  at)  +  -  ^ 

i=0  x0=±a0 


I  I 

XL-l  =  ±a-L-l 


sinh(urx) 


L—l 


=  J” |  cosh(itj  ■  cq) 


i= 0 


Therefore  the  average  over  sinh  functions  must  equal  zero. 

Corollary  1-3: 

z2+”+Z  X  Z  2  X  cosh(uTx)  cosh(vTy) 

Vm- i=±l  *o=±ao  *L-i=±az.-i 

^Z+'  +  sZsZ  -5  Z  cosh(“rj  +  i!7'^ 


yo±i 


yo±i 


yM  =  ±l  ^o  =  ±aO  ^L-1=±«L- 


Proof:  The  proof  is  derived  from  Theorem  1.  It  requires  decomposing  the  cosh 
function  in  a  sum  of  exponentials. 


cosh(ttrx)  cosh(vry) 


_  1  r gUT ic+vT y  _|_  g uTx-vTy  _|_  g -uTx+vTy  _|_  g-uT  x-vT 
4  ^ 


Averaging  over  e±uTx±t7r^  does  not  depend  on  the  sign  of  u  or  v.  Therefore, 
averaging  over  a  product  of  cosine  functions  produces  four  terms  which  averages 
are  equal. 


ZZ 


gUTx+vTy  _ 


ZZ 


e±uTx±vTy 
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^^^^cosh(MTZ)cosh(vTy) 


euTx+vTy 


y  x 


Applying  Theorem  1  and  Corollary  1  to  the  left  side  of  the  expression  proves  the 
Corollary  3. 


z2+  +Z  2  \  2  2  cosh(uTx)cosh(vTy) 

yM=±l  x0=±a0  xL-1=±aL-1 

=  \'L+"*\YJ  5  Z  Z  cosh(“rj  +  i!7'^ 


Xo±l 


yo+i 


yM=±l  *o=±a<)  xL~1=±aL~1 
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APPENDIX  F.  Exact  Likelihood 
Sample  Code  for  Simplification 
of  a  3x3  Hypothesis 

(*  DERIVE  THE  LIKELIHOOD  FUNCTION  *) 

ClearAll[c,  r,  b] ;  Clear [CC,  X,  X ,  njf 

(*  1 .  DEFINE  r  Code  Length  and  Humber  el  Users  +) 

L  =  3;  U  =  3; 

O  z.  construct  t  CDMA  spreading  Matrix  *) 

CC  =  Array [Subscript [c,  &,  tL,  U}]; 

cc  //  MatrlxForm 

'  Co,|}  CijU  C  1,0  1 
Co,l  Ci,!  Ci,i 

^0,2  ^1,2  Cl,  j  , 

(*  3.  CONSTRUCT:  Information  Vector  *) 

B  =  Array [Subscript Lb,  a  -  1,  K]  & ,  (U,  1}]; 

B  //  MatrlxForm 

'  t>6,K  ' 
t>i,K 
1  bi,K  / 

(*  4.  CONSTRUCT:  Correlator  Output  Formatted  as  a  Vector  * ) 

R  =  Array  [subscript  [r,  K,  rtl-1]  &,  {L,  1}]; 

R  //  MatrlxForm 

r*,i 

L  ^,1  j 

(*  5.  COMPUTE:  Correlator  Term  of  the  Conditional  Likelihood  * ) 
x  =  rt,cc.b  t 

X  =  X  [[1,  1]  J 

W,*  (C6,&  rK,0  +  C6,i  rK,l  +  C0,l  rx,a)  + 
bl,K  1  Cl,0  rK,0  +  Cl,l  rK,l  +  Cl,2  rK,i)  4  b2,a  l.C2,0  rK,0  +  C2,l  r(i,l  +  C2,2  rti,2) 
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3  3  3  3 


O  6.  construct i  unnornia I lz cd  signal  model  *) 

S  =  (CC.B); 

S  //  MatriKFom 

b&lK  c&,&  +  bilK  ci|(j  +  bljK  cajfl  ■ 
b&|H  c&,i  +  bi|K  cljL  +  cZxL 
L  c&.a  +  bi|K  ci(  a  +  b2iK  cj, a  , 

O  7,  construct i  Energy  Tern  of  the  Likelihood  Function  *) 
=  (ST.SJ ; 

=  Expand [E^T]  ; 

=  ET  /.  v_x  2  ->  1; 

=  ET[[lr  HI 

®  +  ^  b&|K  biK  c&i(j  cLj&  +  2  beiK  bliK  c0ii  clj  l  +  2  b0fK  biiK  c0il  ciia  + 

^  ^0,K  ^i,K  C0,0  C2f(S  +  ^  ^l,K  ^i,K  C1,<J  C2,0  +  ^  ^0,K  ^2,K  C<S,  1  C2, 1  + 

2  b1(K  b2iK  C1A  e2(1  4  I  b0Hl(  b2(K  c0h2  c2i2  +  2  blH[t  b2<K  c1(2  c2h2 

O  8.  CONSTRUCT:  TOTAL  SQUARED  CORRELATION  #) 

C2  =  ccT.cc; 

V2  =  0; 

For  L  11  s  1,  il  s  B,  il++. 

For  [  12  =  1,  12  s  U,  12++, 
if  [  il  st  12, 

V2  =  V2  +  C2[[il,  ±2]  J  *2 

] 

J; 

]f 

wl  =  Expand [V2]  /,  {_*2  -+  1}  /.  {_'4  -+  1}; 

FML  =  Wl 


(*  9.  CONSTRUCT:  Weight  of  the  Joint  Distribution  of  the  Code 
coefficients  *) 

z  =  Flatten[cc]f 
If  [L  *  1, 

Wl  =  Exp  [  -,0  *  Wl]  f 
For  [  1  =  1,  1  £  Length [z]  r  i++r 
If[!  NumberQ [2 [[!]]], 

rl  =  {  Z  [  [ i]  ]  -+  z]f 
wi  =  wl  i ,  rl; 

Wl  =  5umlWlr  {z,  {-1,  1}]]  f 
]; 


] , 

Wl  =  B*p[-0]  f 
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wi 

LS  +  4  C  ^  g,  C  ^  CSi|  &  CLjL  +  4  Cfljl  ■Gi|[jCSi|5L  +  4  C&jL  CBja  CljL  Ci|2  + 

^  C0H0  C0H1  G2H0  G2h1  4  ^  Gl,0  G1h1  C2,G  G2,l  +  ^  G0,0  C0,2  G2,0  G2H2  + 

^  G1h0  C1h2  G2h0  C2h2  4  ^  GOh1  G0h2  G2h1  G2,2  +  ^  Cl,l  Cl,2  C2,l  G2H2 

32  e  &i'G  +  28 B  e  22$  +  L92  e  s* 


(*  10. 

construct 2  Conditional  likelihood  given  code  matrix  and 
information  vector  *) 

CL  =  Exp  [Sqrt  [2  y]  w  X  -  y  *  ET  -  *  FML]  ; 

£*  11.  COMPUTE t  Average  Likelihood  Function  *) 

Clear [z] ; 

AL  =  CL; 

z  =  Join [Flatten [cc] ,  Flatten[&]]; 

For  [  i  =  1  ,  i  s  Length [Z]  ,  i++. 

If  [  I  Kumbcrp [  Z  [  [i]  ]  ) , 

rl  -  {  Z(  [i]  ] 

AL  =  AL  /  .  rlf 

if  [  Number^ [z  ((!]]  /-  {<=_,_  -0}]p 

AL  =  Slim  [AL,  {  Z  ,  {-I,  1} }  J  , 

AL  =  1  /  2  *  sum  [AL,  {  2  ,  {-I,  1  >  >  ] 

b 

b 

b 

£*  11. 

SIMPLIFY  2  Rules  to  Obtain  the  Likelihood  Function  in  a  Simple  Form  *} 

O  Trivial  Expansion  of  the  Exponents  *) 

ALE  =  AL  /  .  Exp[aaa_]  Exp  [  Expand  [  ZLZkii  ]]  ; 

O  RULEt  Rename  Exp[  r(i]  a  ]  to  P*a  *) 

rule 2  = 

{  Exp  [  Plusf  sMrl _ J  +ejr2 _  a  rR  3 

Exp [Plus fexi] ]*  Subscript [PF  s]  "Times [ex2]  }; 

For  [  i  =  lf  i  s  L ,  i  ++, 

ALS  =  ALS  /.  rule2; 
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O  RULEt  convert  the  variables  PAa  to  Cosh [r (i) *a  ]  +  Sinh[r(i}*a  ]  *) 

rilled  = 

{  Subscript  [P,  s_J  ‘Times  [^2 _ )  1-+ 

Cosh  [Times  [ex2]  *  Subscript [  r r  KP  s]]  + 

Sinb[Times  [ex2]  *  Subscript  [r,  Kr  is)  )  }) 

ALS  =  Expand [AL5  /,  rule!]? 

O  12 .  FINAL  LIKELIHOOD  FORM  *) 

als  =  Expand  [als /wi  ] 


24  e  54  ^  Jt  Cosh  [V  2  V7  rKli}]  Cosh  [V  2  V7  rKll]  Cosh  [V  2  VT  rKjl] 


144  e  22  ^  3  r  Cosh 

32  e  54^  +  288  e  S2- 

'[VI  V7  ^Kjd]  Cosh 

3  4-  192  £  Gl'3 

j  V  2  V7  rKji]  Cosh  [ 

-  + 

V  2  V7  rKj2] 

48®  G*'3TCosb[-\ 

32®  54*  +  2 8 8  ®  2 
/  2  Vy"  rxjO  ]  Cosh  [V 

:2*  +  192  e  G* 

2  Vy  rK(i]  Cosh 

1 

v'2  Vy"  rKja] 

24®  22 &  11  y  Cosh 

32®  54*  +  288  ®  22 

[  3  V  2  V7  rKjfl]  Cos] 

^  +  192  e  G^ 

h.  [V  2  V7  rKji]  Cosh 

[VI  V7  rKj3] 

48®  G*  llTCosh[ 

32®  541,5  +  283  e 

3  V  2  Vy  rKj0]  Cosh 

22l3+  192®  Gl? 

[VI  V7  rKji]  Cosh  [ 

VI  V7iKiS] 

24®  22(i  llT  Cosh 

32®  54*  +  288  ®  22^  4  192  e 
[V  2  Vy  rKjfll  Cosh  [3  VI  Vy  rKli]  cosh 

[VI  V7  rK(3] 

48®  G^  llT  Cosh 

32®  54l,J  4  283  ®  22*  4-  192  e  ** 

V  2  Vy"  J^k.4  ]  Cosh  [  3  V  2  Vy  1  ]  Cosh  [ 

VI  VyiK(J] 

24®  22  &  19  T  Cosh 

32®  54*  +  288  e  2i*  4  192  e  Gl'5 
[  3  V 2  V7  ]  Cosh  [  3  V  2  Vy"  1  ]  1 

Cosh 

[VI  Vy  i'Kji  ] 

24®  22!i  11 T  Cosh 

32®  G4e  +  288  ® 

[v  2  V7  *Kti\]  Cosh  [ 

+  192  «  4 ^ 

V  2  Vy  rKji  ]  Cosh  [3  V  2  V7  rKji] 

48®  G^  llr  Cosh  [ 

32®  541,5  4  288  ® 

V  2  V7  **,0]  Cosh  ['■ 

22l3+  192®  G* 

/  2  V7  rKj  1  ]  Cosh  [  3  V 2  V7  *Kji  ] 

24®  22/1  19  T  Cosh 

32  e  54*  +  288  e  : 

[3  VI  V7  \  Cos] 

!2^  +  192  ®  G* 

ti  [V  2  Vy  1  ]  cosh 

[3  VI  Vy"  ^2] 

24®  22  19  T  Cosh 

32®  54^  +  288  ®  3i*  4  192  e  Gl9 

v  2  Vy  Cosh  [3  VI  Vy"  rKji]  Cosh 

[3  VI  Vy"  ^i] 

32®  54e  +  288  e 
8®  si r  Cosh  [ 3  VI"  V7  rKjll]  Cosh 

22*  +  192®  ** 

[  3  V 2  V7  rKj  1  ]  Cosh  [  3  V 2  V7  rKja  ] 

32  e  5Jj3  +  2B8  e  22^  +  192  e 
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(*  13.  case:  when  $=y  *) 

LF  =  ALS  /,  {/3  y} 

24  6  51  ^  Cosh[V  2  Vy  rKi&]  Cosh  [V  2  Vy  rKjL]  Cosh  [V  2  Vy  rKi2] 
-  + 

32e  “*  +  238  e  22  *  +  192  e 

144  e  25  v  Cosh  [V  2  Vy"  rKjoj  Cosh [VI  Vy"  rKjJ  |  Cosh  [v  2  Vy"  rKj2] 
32  e  5Jlf  +  2B8  e  23  T  +  192  e  fiT 


48*  9  T  Cosh 

V  2  Vy  rKjD]  Cosh  [ 

V  2  Vy"  ^Kji]  Cosh 

[v'2  VT  rKjj] 

32  e  *  2B3  e  22  *  +  192  e 

24*  33  y  Cosh  [3  VT  Vy"  rKjfl]  Cosh[VT  Vy"  rKji]  Cosh 

v  2  Vy  rKja  ] 

32  e  5Jt  +  2B8  e 

48*  17 y  Cosh  3V2  Vy  Cosh 

22t-  192  e  filf 

[VI  Vy  rKjl]  Cosh 

[V 2  Vv  rKj2] 

32  e  +  2B8 

24  *  33 Cosh  [VI"  V7  rKj(i]  Cosh 

e  22 T+  192  e 

[  3  V  2  Vy  rKji]  Cosh 

v  2  Vy  rKji] 

32  *  *  2B8 

48*  17  y  Cosh  V  2  Vy"  rKjo|  Cosh 

*  22*  +  192  e  6r 
[3  v'2  Vy  rKii]  Cosh 

v'2  Vy  rKj2] 

32  e5i!f  +  2BSe 
2  4*  91  ■'  Cosh  [  3  V 2  V7  rKj0  ]  Cosh 

22*-  192  * 

[  3  V 2  Vy"  rKji]  Cosh 

[VI  Vy  rKja] 

32  *  5Jt  +  2Bi 

24  *  33  y  Cosh  V  2  V7  rKj  a  |  Cosh 

i  e  22  *  +  192  e 

[V  2  Vy"  rK,i]  Cosh  [3  v  2  Vy  rKia] 

32  e  2B8 

48*  17  T  Cosh  V  2  V7  ^KjD  Cosh 

e  22 T  +  192  e  6f 

[  V  2  V7  rKj  1  ]  Cosh  [  3  V  2  Vy  ^2  ] 

32  e  54*  *  2B8  e 
24*  41 T  Cosh  [3  VT  V7  rKjll]  Cosh 

23  *  -  192  *  6* 

[VI  VT  rKji]  Cosh 

I 

[3  v'2  VT  rKji  ] 

32  *  5Jt  +  268  e  22t  +  192  e 

24*  41 T  Cosh  V 2  V7  rKjo|  Cosh  [3V  2  Vy"  rKjl]  Cosh 

[3  v'2  V7  rKi2] 

32  e  54^+  2B3  e  22^  +  192  e  ^ 


8e  11  T  Cosh  [  3  V  2  W  ]  Cosh  [  3  V  2  Vy  rKj  i  ]  Cosh  [  3  V  2  Vy  rKji  ] 

32*  slT  +  288  e  22t  +  192  *  ‘T 
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(*  14.  CONSTRUCT:  EMPIRICAL  LIKELIHOOD  *) 


(»  Multi combi nations  *) 

BZ  =  Floor  [V/2]  +1; 

MDigitBascBZofN  =  Function  [{W,  BZf  Af } ,  1st  =  { 0}  f 
if  [  BZ  >  1,  1st  =  integer  Digits  [W,  BZ]]} 

if  l  Af  >  Length [ 1st]  ||  Length [  1st  ]=  0,  0,  1st  [  [Length [1st]  +  1  -  AT]  ]]  ] 


aiL-l 

A=  £  *[l,U, 

H-& 


1 

L 


2 


( (L  +  1)  4 (MDigltBaseBZofN [N, 


BZ  ,  M2])  -  1 )  | 


Li 

|  |  lExp  l-y  (2  + MDlgitBa&CBZOf N  [N,  BZ ,  M]  +  Mod  [U,  2]  )  '2/1] 

Cosh  (2  * MDigltBaseBZofN [N,  BZ,  M]  +Mod[U,  2])  y2  y  /  I  rKi^ 


Mil 


e  J>  Cosh  [ VT  V  y  rKi&]  Cosh  [V  2  V7  rKiL]  Cosh  [V  2  V7  rKil]  a[3,  3,  0]  + 
e  11  *  Cosh  [3  V  2  W  rK,o  [  Cosh  [V  2  Vy"  rK(1  |  Cosh  V  2  Vy"  rK(2]arir  3r  1]  + 
e  11  *  Cosh[vT  V7  rK(fl]  Cosh  [3  VT  V7  rK#1]  Cosh[vT  V7  rK(2]  a  [3,  3P  11  + 
e  11  v  Cosh  [  V  2  Vy"  rK,o  II  Cosh  V  2  Vy"  iK,i  i  Cosh [3 V 2  Vy"  xK#2  ]  a[3P  3P  1]  + 
e  11  v  Cosh  [3  V2"  V7  rK(il]  Cosh  [3  V2"  V7  rK#1]  CoshfVT  Vy  r^2]  a[3P  3,  2)  + 
e  v  Cosh  [  3  V  2  Vy"  *k,o  I  Cosh  [V  2  Vy"  rK(  j  ]  Cosh  [  3  V  2  Vy"  rK(2  ]  a  f  3 ,  3  p  2 1  + 
e  1S  *  Cosh  [V"2  W  rK(0  ]  Cosh  [3  V2"  W  rK(1  |  Cosh  j  3  V  2  Vy"  rK(2]of3P  3,  2)  + 
e  2T*  Cosh  [3  VI"  V7  r^a]  CoshflVT  V7  rK#1]  Cosh[3V"2  V7  rK#a]  a[3P  3r  3] 


(*  15.  CASE:  When  for  Different  Feature  Vectors  *) 

*[3,  3,  0]  * 

(als  //.  |cosh[VT  V7  rK,o]  Cosh [v 2  Vy"  rK,i]  cosh [V 2  Vy  r^jj  -+  y, 

Cosh  [ _ ]  -»  0  |  J  /  .  {y  -j  1 ,  y  ^  Q} 

Plot [Log [n [3,  3,  0]],  {fl,  0,  20},  PlotRangc  -»■  Full, 

PlotLabel  -+  [0]  ”,  AKesLabel  -+  {"0”  P  "Log  [a3i3f0  [jSJ  J  ”  }] 

Together  [  n[3,  3f  0]  /.  {  Exp  [  u _ z  *  (-m/  16)  }  ] 

24  e  MB 

32  e ■  +  2BB  e  32  ^  +  192  e 

144eS2S  4Be6'3 

32  e  54^  4  288  e  22fl  4  192  e  **  32  e  +  266  t  12 ^  +  192  e 
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-1.0  ' 


-I  I  - 


-1.2  ■ 


L0 


L5 


20 


3  (2  +  62  +  zJ) 

4  (6  4-  9  2  +  z1) 
n[3,  3,  1]  = 

^als  // .  |  cosh  1 3  Vi"  Vr  tK,o ]  cosh  |  Vi"  Vr  r*r  t  J  cosh  | Vi"  Vr  rK,2 ]  -+  T r 

coshf _ ]  -*  o  })  / .  {y  -*  1,  t  -*  o> 

Plot [Log [n [3,  3,  1]J,  in,  0,  20},  PlotRange  -»■  Full, 

FlotLabel -►  AKesLabel -►  {’'jS"'  F  "'Log  1  Ifii  J ] 

Together  [  a  [3,  3,  1]  / .  {  Exp  [ti _ ft]  -+  z  *  (-11/  16)  }  ] 

24  e  12fh  48efifl 

32  e  +  2BB  e  22 -G  +  L 92  e  t-G  +  3 2  e  *  +  2 83  e  22  6  +  1 91  e  fij3 

Ojjj  [J5J 

I^II 

-1.4  :f" 


-1.5 


-l.fi 


-1.7 


-L.S 


—i  1_ 1_ 1  L 


m 


[5 


20 
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*[3,  3,  2]  = 

Together | 

|als  f  f  -  {  cosli^3  V2"  Vt”  rK0  cosh^3  V2"  Vr"  rn,i]  coshj^V^"  W  rKrz]  -*  Vr 

Cosh  [ _ -+  0  /  .  {.JT  -►  1 ,  y  -►  0}  J 

Plot  [Log  [a [3,  3,  2]  ]  ,  {fl,  0,  20},  Plot Range  -+  Full, 

P  lot  Lab  el  -+  "<*3,3,2  [£]  "  ,  AxcsLabel  -+  {"£”  r  "Log  [<*3,3,2  [£]  J  "}  ] 

Together  [  n[3,  3,  2]  / .  {  Exp  [  11 _ £]  -+  z  ‘  (-u  /  16) }  ] 

]eJ1^ 

4/1  +  9  e12*3  +  6  eja^  ) 


3  z 

4  f  6  +  9  z  +  zJ) 

n[3,  3,  3]  = 

Together | 

|als  / / .  |  Cosh^3  V2"  V7  rn,o]  Coshj  3  Vi"  V7  rnri]  Coshj^S  -^[T  W  rK  jJ  ■*  Yj 

Cosh  [ _ J  -+  0  })  /.  u+lfT+B}] 

Plot [Log [n [3,  3,  3]],  {0,  0,  20},  PlotRangc  -+  Full, 

PlotLabel  -+  ”<*3,3,3  [£]” ,  AxcsLabel-*  "Log  [<*3,3,3  [^3 J  | "}  ] 

Together  l  n[3,  3,  3]  /.  {  Exp  [u _ fl]  -+  z  '  u  /  16 )  }  ] 


# 


Created  with  Wolfram  Mathematica  7.0 


Approved  for  Public  Release;  Distribution  Unlimited 

43 


APPENDIX  G.  Calculation  of  aL :U,ci(P) 

Using  Programming  Code 


The  implementation  of  the  average  likelihood  takes  the  following  form 

L— 1 

A(r|Jf)  =  ^  aL,u,aY\e~r^=ola^c  osh(72y  ■  Re{rt}  ■  a£) 

All  a  i= 0 


where 


aL,U,a 


2Sum(a^=0)  'SjAU  r  Sce 


Q  ^’T(C) 


2j/U/  a  2Sum(a;=0)  ^-iAll  T  Sce  ®  P  (  ) 


A  Matlab  implementation  of  the  parameters  is  shown  below.  The  code  was  executed 
using  AFRL  Condor  HPC  cluster  for  calculating  the  values  of  alphas  when  /?  -»  oo. 


%  Calculation  of  the  numerator  of  the  alpha  parameters  using  formula. 

%  A.  Vega,  Last  Modified:  October  8,  2013 
%  inputs:  a  =  vector  vector,  U  =  number  of  users 
%  output:  matrix  =  [  TSC,  occurrences  ] 

%  formatted  output  =  occurrences  *  eA  (beta*TSC) 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

function  output  =  alphaParam (a, U) 


%  Ensure  positive  coefficients 
a  =  abs  (a) ; 


%  Code  length  of  the  spreading  matrix 
L  =  length (s) ; 

%  Calculate  possible  column  permutations  such  that  sum (cols)  >=  0 
col Perm  =  []; 

colPerm ( 1 )  . x  ( 1 ,  : )  =  ones(l,U); 
for  kl  =  l:floor(U/2) 
x  =  ones ( 1 , U) ; 
k2  =  kl; 
while  (  k2  >  0  ) 

x(l,k2)  =  -1; 
k2  =  k2  -  1  ; 

end 

colPerm ( kl+1 ). x  =  unique (perms (x) , ' rows ' ) 

end 


%  Display  possible  row  vectors 
if  (  0  ) 

for  k3  =  1 : size  (a ( : ) , 1 ) 
squeeze (a (k3)  . x  (  : ,  : ) ) 

end 

end 


%  Construct  matrices  using  row  vectors 
group  =  abs (U-a) /2+1 ; 
elems  =  zeros (1,U); 
for  kl  =  1 : L 
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elems(l,kl)  =  size (colPerm (group (kl )). x, 1 ) ; 

end 

%  This  algorithm  is  inefficient  for  L>6 
tic 

stack  =  [];  mtx  =  zeros  (L,U) ; 

start  =  ones(l,L); 

while (  all (start<=elems)  ) 

%  Generate  Matrix 
for  kl  =  1 : L 

mtx(kl,:)  =  colPerm (group (kl) ). x (start (1, kl) , :) ; 

end 

%  Calculate  Total  Square  Correlation 

tau  =  0; 

h2  =  mtx'*mtx; 

for  kl  =  1:L 

for  k2  =  1 : L 

if (  kl  ~=  k2  ) 

tau  =  tau  +  h2(kl,k2)A2; 

end 

end 

end 

%  Store  value 
stack  =  [stack, tau] ; 


%end 

%  This  algorithm  explores  all  possible  combinations  of  a  row  in  a  matrix 
%  The  sum  of  the  row  >=  0 
k2  =  1; 

while (  k2  <=  L  ) 

c  =  start (l,k2)  +  1; 
if (  c  >  elems(l,k2)  ) 
if (  k2  <  L  ) 

start ( lr k2 )  =  1 ; 
k2  =  k2  +  1  ; 

else 

start ( lt  k2 )  =  c; 
k2  =  k2  +  1  ; 

end 

else 

start ( 1 , k2 )  =  c; 
k2  =  L+l; 

end 

end 

end 

toe 

%  Count  all  unique  values 
vals  =  unique (stack' , 'rows' ) ' ; 
svals  =  zeros (size (vals) ) ; 
for  kl  =  1 : length (vals ) 

svals (l,kl)  =  sum (stack==vals (kl ) ) ; 

end 

%  Display  results 

[vals ' ,  svals ’ *2 AL/2 A (sum (s==0 ) ) ] 
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APPENDIX  H.  Various  aL  U s(fi)  for  L  <  4 


The  coefficients  a  depends  on  the  dimensions  L  (rows)  and  U  (columns)  of  the 
spreading  matrix  C.  Note  that  aL  Urf(J3)  =  aL(J  ^(/?)  if  a  is  a  permutation  of  b. 

L  =  1,  U  =  1 

=  1 

L  =  2,  U  =  1 

★  _  ^ 

^2,l,[l,l]r  —  ^ 

L  =  2,  U  =  2 
z  = 

a2,2,[0,0]T  =  Z/(2  *  C1  +  Z)) 
a2,2,[2,0]T  ~  1/(2  *  (1  +  z)) 
a2,2,[2,2]r  =  Z/(2  *  C1  +  Z)) 

L  =  3,  U  =  1 

a3,l,[l,l,l]T  = 

L  =  3,  U  =  2 

Z  —  g— 16- P 

a3,2,[0,0,0]T  =  Z/(2  *  (3  + 
a3,2,[2,0,0]r  —  1/(2  *  (3  +  Z)) 
a3,2,[2,2,0]r  =  1/(2  *  (3  +  Z)) 
a3,3,[2,2,2]T  =  Z/(2  *  (3  +  Z)) 

L  =  3,  U  =  3 

z  =  e~16'P 

a*  2  [i  i  i]T  =  3  *  (2  +  6z  +  z3)/(4  *  (6  +  9z  +  z3)) 
a3,2,[3,i,i]r  =  3  *  (2  +  z)/(4  *  (6  +  9z  +  z3)) 
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a3,3,[3,3,l]r  =  3z/(4  *  (6  +  9z  +  Z3)) 
a3  3  [3  3  3]T  =  Z3/(4  *  (6  +  9z  +  Z3)) 

L  =  4,  U  =  1 

a4,l,[l,l,l,l]T  =  1 
L  =  4,  U  =  2 
z  = 

a4,2,[0,0,0,0]T  =  z4/(2  *  (!  +  z)2  *  (3  -  2z  +  z2)) 
a4,2,[2,0,0,0]T  =  Z/(2  *  C1  +  z)2  *  (3  -  2z  +  z2)) 

a4,2,[2,2,0,0]r  =  1/(2  *  (!  +  z)2  *  (3  -  2z  +  z2)) 

a4,2,[2,2,2,0]r  =  z/(2  *  (1  +  z)2  *  (3  -  2z  +  z2)) 

a4,2,[2,2,2,2]T4  =  ZV(2  *  (l  +  z)2  *  (3  -  2z  +  Z2)) 

L  =  4,  U  =  3 
z  = 

a4 3  [i  1 1  jjt  =  3z  *  (12  +  6z  +  8z2  +  z5) /(4  *  (6  +  36z  +  9z2  +  12z3  +  z6)) 

a* 3  [31u]r  =  3  *  (2  +  6z  4-  z3)/(4  *  (6  +  36z  +  9z2  +  12z3  +  z6)) 

a4  3  [3  3 1  i]T  =  3z  *  (2  +  z)/(4  *  (6  +  36z  +  9z2  +  12z3  +  z6)) 

a4  3  [3  3  3  ijT  =  3z3/(4  *  (6  +  36z  +  9z2  +  12z3  +  z6)) 

a4  3  [3  3  3  3]T  =  z6/(4  *  (6  +  36z  +  9z2  +  12z3  +  z6)) 

L  =  4,  U  =  4 

z  =  e-W 

q(z )  =  (8  *  (6  +  96z3  +  144z4  +  144z7  +  81z8  +  24z12  +  16z15  +  z24)) 

a4,4,[o,o,o,o]T  =  3z4  *  (12  +  6z4  +  8z8  +  z 20)/q(z) 

a4,4,[2,0,0,0]T  =  12z3  *  (2  +  6z4  +  z12)/q(z) 

a4,4,[2,2,0,0]T  =  12z4  *  (6  +  5z4  +  z8)/q(z) 

a4,4,[2,2,2,0]T  =  12z3  *  (6  +  9z4  +  z12)/(?(z) 
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a4,4,[2,2,2,2]r 

★ 

a4,4,[4,0,0,0]r 

a4,4,[4,2,0,0]T 

a4,4,[4,2,2,0]T 

a4,4,[4,2,2,2]r 

a4,4,[4,4,0,0]r 

a4,4,[4,4,2,0]r 

a4,4,[4,4,2,2]r 

a4,4,[4,4,4,0]r 

a4,4,[4,4,4,2]T 

a4,4,[4,4,4,4]r 


4  *  (6  +  36z4  +  9z8  +  12z12  +  z24)/q(z) 
3  *  (2  +  6z4  +  z12)/q(z ) 

12z3  *  (2  +  z4)/q(z) 

12z4  *  (3  +  z4)/q(z) 

4z3  *  (6  +  9z4  +  z12)/g(z) 

3z4  *  (2  +  z4)/g(z) 

12z7/g(z) 

4z8  *  (3  +  z4)/q(z) 

3  z12/q(z) 

4z15/q(z) 

z24/q(z) 
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APPENDIX  I. 


Limit  of  aL'U,ci(.P)  for 
General  Cases 


We  define  the  parameter  cq,,o,a(/?)  using: 

2l  q  ^ 

2Sv.m(ai=0)  ^-iAll  T  Sze  S-£r  ® 

(Xl.uAP)  = - iZ - : - — 

^iAll  a  2Sum(a;=0)  ^-iAll  T  2jZE  .->gT  ®  ^ 

Proposition  5.  The  limit  of  ar,u,a(/?)  as  /?  approaches  infinity  depends  on  the 
feature  vectors  as  follows: 


lim  ccLiUmS  (jff)  = 

/?-»  oo 


0  /or  a  =£  d* 


2Sum(ai=0 ) 


>  0  /or  a  =  a* 


uAll  a  2Sum(ai=o) 


Proof:  The  expression  can  be  converted  in  a  polynomial  rational  function  by 
substituting  e~P  -*  y  and  divide  by  the  exponent  with  the  minimum  Total  Square 
Correlation  yTmin  jn  both  side  of  the  equation. 

2l 

2Sum(at=0)  ^jAU  t  2ze  S-^t  V 

«L,u,a(P)  = - 2t - - - 

^ iAll  a  2 Sum(af=o )  ^ iAll  T  5jZE  S-^T  V 
2l 

2Sum(aj=o)  ^ iAll  T  5jZE  S-^T  y  m^n 

~  2^ 

T.AU  a  2Sum(.ai=0)  ^All  T  2ze  SUt  yT~Tmln 

T  =  T (Z)  >  Tm;n  >  0 
a  =  abs(Z  ■  1) 

Since  r  —  Tmin  >  0,  taking  the  limit: 

t-t  ■  f0  for  x  =£  Tmin 
lim  yT  Tmin  =  ]  ' 

y^O-7  (1  /or  T  =  Tmin 

eliminates  all  terms  with  Total  Square  Correlation  terms  greater  than  Tmin.  The 
result  is: 


aL,U,a(/0  — 


2sum(ai= 0)  ^zesn 


uAll  a  2^um{ai=o) 


ZzE^ 


>  0 
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The  summation  over  all  the  terms  in  the  set  is  just  the  cardinality  of  the  set, 
therefore: 


a 


l,u , 


= 


2Sum(ai=0 ) 


uAll  a 


2Sum(ai=0 ) 


CL,Tyi 
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APPENDIX  J.  Feature  Vectors 
Code  length  L  <  6 


Table  J-l.  Several  Feature  Vectors 


L=2 

L=3 

L=4 

L=5 

L=6 

U=1 

[1,1]T 

[ixiy 

[1,1,1,1]T 

[1,1,1,1,1]T 

[1,1,1,1,1,1V 

U=2 

[2,0]^ 

[2,0,0V 

[2,2,0V 

[2,2,0,0]T 

[2,2,0,0,0]T 

[2,2,2,0,0]T 

[2,2,2,0,0,0V 

U=3 

[3,1,1V 

[1,1,1V 

[3,3,1,1]T 

[3.3.1.1.1] T 

[3.3.3.1.1] T 

[3,1,1,1,1,1V 

[3,3,1,1,1,1V 

U=4 

[2,2,2,2]T 

[4,0,0,0]T 

[2,2,2,2,2]T 

[4,0,0,0,0]T 

[2,2,2,2,0]T 

[4,2,0,0,0]T 

[4,2,2,2,2]T 

[4,4,0,0,0]T 

not 

calculated 

U=5 

[3,3,3,3,3]T 

[5.1.1.1.1] T 

[3.3.1.1.1] T 

not 

calculated 

U=6 

[2,2,2,0,0,0]T 

[4,2,2,2,0,0]T 

The  Corresponding  Lowest  TSC  Spreading  Matrices 

Case:  L  =  2,  U  =  2,Tmin  =  0 


C*  = 


c 


3 


Case:  L  =  3,  U  —  2 ,Tmin  —  2 
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c* 


■+ 

+' 

■+ 

+■ 

+ 

— 

f 

+ 

— 

.+ 

+. 

+. 

Case:  L  = 

3  ,U 

= 

3, 

^min 

= 

6 

'+ 

+ 

+' 

'+ 

+ 

C*  = 

+ 

— 

+ 

t 

+ 

— 

+ 

.+ 

+ 

+ 

+. 

Case:  L  =  4,  U  =  2,  Tmin  =  0 


C* 


-+  +' 

+  - 

+  + 

-+  -- 


Case:  L  =  4,  U  =  3,  Tmin  =  0 


C* 


[+  +  +] 
+  -  + 
-  +  + 

L+  +  -J 


Case:  L  =  4,  U  =  4,  rmin  =  0 


-+  +  +  +- 

-+  +  -+- 

+  -  +  - 

+  -  +  + 

+  +  -  - 

-  +  +  + 

L+  +  -  +J 

-+  +  +  -- 

Case:  L  =  5,  U  =  2,  r min  =  2 


C* 


-+  +- 
+  - 
+  + 
+  - 
-+  +- 
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Case:  L  =  5,  U  =  3,  r min  =  6 


r  +  +i 

r+  +  +i 

+  +  + 

+  +  - 

+  -  + 

p 

-  -  + 

+  +  - 

+  +  + 

1+  -  +J 

L+  +  +J 

Case:  L  =  5,  U  =  4,  rmi„  =  12 


r+  + 

+  +i 

r+  +  -  +i 

r+  +  -  +i 

- 

+  -  +  - 

+  +  -  - 

t 

+  -  +  + 
-  +  +  + 

/ 

+  -  +  + 
-  +  +  + 

/ 

+  - 

-  + 

+  +  +  + 

+  +  +  + 

-+  + 

+  +. 

L+  +  +  +J 

L+  +  +  +J 

- 

C*  = 


+  +  +  + 

+  +  +  + 
+  +  +  + 

+  +  +  + 
+  +  +  + 


+  +  -  + 

+  -  +  + 
-  +  +  + 

+  +  +  - 
+  +  +  + 


Case:  L  =  5,  U  =  S,Tmin  =  20 


-+ 

— 

+ 

+ 

+- 

-+ 

+ 

— 

+ 

+- 

~+ 

+ 

+ 

+ 

+- 

— 

+ 

+ 

+ 

+ 

+ 

— 

+ 

+ 

+ 

— 

+ 

+ 

— 

+ 

+ 

+ 

— 

+ 

— 

t 

— 

+ 

+ 

+ 

+ 

) 

— 

+ 

— 

+ 

+ 

+ 

+ 

+ 

— 

— 

+ 

+ 

+ 

+ 

+ 

— 

— 

+ 

+ 

+ 

-+ 

+ 

— 

— 

+- 

-+ 

+ 

+ 

+ 

— - 

+ 

+ 

+ 

— - 
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APPENDIX  K.  Simplification  Using 
Sylvester's  Construction 


The  average  likelihood  function  under  the  assumption  of  a  Sylvester’s  construction 
is: 


A(r\M  =  Hl)  = 


~(r,r) 


L—l 


SnNo)L 


2  Y\e  ylfi|2  cosh(V2y  -n). 


All  permutations  i= o 
of  f  *  elements 


The  likelihood  ratio  between  hypotheses  Jf2L  =  H2 l  and  =  HL  will  be 

considered: 


A(r\K  =  H2L )  = 


,-(r,r) 


-  2L-1 

i  2 


SnNo) 


2  L 


^  e  Y l2Ll2  cosh(-y2y  ■  2 L  •  r {) 


k= 0 


and 


2(f  \H  =  Hl)  = 


SnNo)L 


L—l 

I' 

k= 0 


-Y  M2 


cosh(A/2 y  ■  L  -ri) 


For  each  sample  vector  r2Lxl  in  the  H2L  there  are  two  sample  vectors  rLxl  in  KL. 
These  two  vectors  are  statistically  independent,  so  the  likelihood  in  the 
denominator  is  multiplied. 


A (r0,ri,  ...,r2L) 


_ 2(r0,r1,...,r2L|Jf  =  tf2L) _ 

2(r(),  r:, ... ,  =  //L)  ■  2(ri; ri+1, ... ,  r2L \ J-C  =  HL) 


Substituding  in  the  previous  equation  gives: 

A(r0,ri,  ...,r2L) 

((^^)2 '  Zk=l  e“K|2L|2  cosh(y2y  ■  21  ■  a) 


(t^¥)2  '  e  KlL|2  cosh(V2y  ■  l  ■  A)  ■  ( ■  2k= o  e  y|L|2  cosh(72^  ■  L  ■  ri+L) 


re-(r2>r2)\ 2 
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Note  that  (A,  A)  and  ( r2,r2 )  is  the  dot  product  between  the  first  and  second  sample 
of  0fL.  The  sum  of  both  dot  products  equals  the  dot  product  in  the  numerator. 

(A>  A)  +  {A- A)  =  <r,r) 

So,  the  product  between 

11  1 
/e~(riA)\2  /e-<r2,r2)\2  /  e~(f,f)  \2 

cancels  with  the  same  term  in  the  numerator. 

The  final  form  is  given  by: 

_ _ Sfcio1  cosh(72y  ■  2 L  ■  Re{rk}) _ 

KA^l)  YJf=  o  cosh(72y  ■  L  ■  flefa})  ■  1^7^  cosh(72y  ■  L  ■  fle{rfc}) 


Using  the  double  angle  formulas  for  cosh  functions  helps  reducing  the  amount  of 
calculations  by  reusing  the  terms  in  the  denominator. 

ACgjgzL)  giogdf  -  1) 

A(r|WL)  Ifc=0  dk  ■  YAAl+i  dk 


dk  =  cosh  (^27  1  L  '  Mrfc}) 


Note:  It  was  discovered  that  the  likelihood  ratio  is  only  valid  when  the  energy  of  the 
signal  is  assumed  to  be  known.  A  more  correct  expression  will  consider  the  effect  of 
the  dependency  of  the  hypothesis  on  the  energy  of  the  CDMA  signal.  In  other  words, 
dk  is  a  function  of  Of. 
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List  of  Acronyms 


ALRT 

CDMA 

LRT 

MPI 

Re{} 

TSC,  r 
(8) 

C*,i 

y 

Vs 

A  2L<7) 

E 

A(r\K) 
A(r\K,  C,  b ) 
C 
C * 
b 

aL,U,a 

a* 

a* 


Average  Likelihood  Ratio  Test 
Code  Division  Multiple  Access 
Likelihood  Ratio  Test 
Multiple  Programming  Interface 
Real  part 

Total  Square  Correlation 

Kronecker  product 

Column  vector 

Chip  Signal  to  Noise  Ratio 

Symbol  Signal  to  Noise  Ratio 

Likelihood  ratio  [Hypothesis  2L  versus  L) 

Average  energy  per  symbol 

Likelihood  function 

Conditional  likelihood  function 

Spreading  matrix 

TSC  optimal  spreading  matrix 

Information  vector 

Coefficients  of  the  CDMA  likelihood  function 
Non-vanishing  alpha  coefficients 
Feature  vector  for  CDMA  detection 
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